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MULTIVARIATE  NONPARAMETRIC  CLASSES 
IN  RELIABILITY 


by 

Henry  W.  Block  and  Thomas  H.  Savlts 


ABSTRACT 

Recent  research  in  multivariate  nonparametric  classes  in  reliability 


is  surveyed 


1.  Introduction.  This  paper  is  a  sequel  to  the  survey  paper  of  Hollander 


and  Proschan  (1984)  who  examine  univariate  nonpar ame trie  classes  and  methods 
in  reliability.  In  this  paper *ve  will  examine  multivariate  nonparametric 
classes  and  methods  in  reliability. 

Hollander  and  Proschan  (1984)  describe^the  various  univariate  nonpara¬ 
metric  classes  in  reliability.  The  classes  of  adverse  aging  described  include 
the  IFR,  I FRA,  NBU,  NBUE  and  DMRL  classes.  The  dual  classes  of  beneficial  aging 
are  also  covered.  Several  new  univariate  classes  have  been  introduced  since 
that  time.  One  that  we  will  briefly  mention  fis  the  HNBUE  class,  since  we  are 
aware  of  several  multivariate  generalizations  of  this  class. 

The  univariate  classes  in  reliability  are  Important  in  applications  con¬ 
cerning  systems  where  the  components  can  be  assumed  to  be  independent.  In 
this  case  the  components  are  often  assumed  to  experience  wearout  or  beneficial 
aging  of  a  similar  type.  For  example,  it  is  often  reasonable  to  assume  that 
components  have  increasing  failure  rate  (IFR).  In  making  this  IFR  assumption 
it  is  implicit  that  each  component  separately  experiences  wear  and  no  inter¬ 
actions  among  components  can  occur.  However  in  many  realistic  situations,  adverse 
wear  on  one  component  will  promulgate  adverse  wear  on  other  components.  From 
another  point  of  view  a  common  environment  will  cause  components  to  behave  sim¬ 
ilarly.  In  either  situation,  it  is  clear  that  an  assumption  of  independence 
on  the  components  would  not  be  valid.  Consequently  multivariate  concepts  of 
adverse  or  beneficial  aging  are  required. 

C 

Multivariate  nonparametric  classes  have  been  proposed  as  early  as  1970. 

For  background  and  references  as  well  as  some  discussion  of  univariate  classes 

with  multivariate  generalizations  in  mind  see  Block  and  Savits  (1981).  In 

the  present  paper  we  shall  only  describe  a  few  fundamental  developments  prior  to  1981 
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and  focus  on  developments  since  then.  The  coverage  will  not  be  exhaustive  but 
will  emphasize  the  topics  which  we  feel  are  most  important. 

Section  2  deals  with  multivariate  nonparametric  classes.  In  section  2.1 
multivariate  IFRA  is  discussed  with  emphasis  on  the  Block  and  Savits  (1980) 
class.  Multivariate  NBU  is  covered  in  Section  2.2  and  multivariate  NBUE 
classes  are  mentioned  in  Section  2.3.  New  developments  in  multivariate  IFR  are 
considered  in  Section  2.4  and  in  Section  2.5  the  topics  of  multivariate  DMRL 
and  HNBUE  are  touched  on. 

Familiarity  with  the  univariate  classes  is  assumed.  The  basic  reference 
for  the  IFR,  IFRA,  NBU  and  NBUE  classes  is  Barlow  and  Proachan  (1981).  See 
also  Block  and  Savits  (1981).  For  information  on  the  DMRL  class  see  Hollander 
and  Proachan  (1984).  The  HNBUE  class  is  relatively  recent  and  the  best  ref¬ 
erences  are  the  original  articles.  See  for  example,  KlejsjB  (1982)  and  the 
references  contained  there. 

2.  Multivariate  Nonparametric  Classes.  Many  multivariate  versions  of  the 
univariate  classes  were  proposed  using  generalizations  of  various  failure  rate 
functions.  These  multivariate  classes  were  extensively  discussed  in  Block  and 
Savits  (1981).  Other  classes  were  proposed  by  attempting  to  imitate  univariate 
definitions  in  a  multivariate  setting.  (See  also  Block  and  Savits  (1981).)  One 
of  the  most  Important  of  these  extensions  was  due  to  Block  and  Savits  (1980) 
who  generalized  the  IFRA  class.  This  multivariate  class  was  proposed  to  parallel  the 
developments  of  the  univariate  case  where  the  IFRA  class  possessed  many  important  closure 
properties.  As  in  the  univariate  case  the  following  multivariate  class  of  IFRA, 
designated  the  MIFRA  class,  satisfies  important  closure  properties.  First,  as  in 
the  univariate  case,  monotone  systems  with  MIFRA  lifetimes  have  MIFRA  lifetimes 
and  independent  sums  of  MIFRA  lifetimes  are  MIFRA.  From  the  multivariate  point 
of  view,  subfamilies  of  MIFRA  are  MIFRA,  conjunctions  of  Independent  MIFRA  are 
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MIFRA,  scaled  MIFRA  lifetimes  are  MIFRA,  and  various  other  properties  are  satisfied. 

We  discuss  this  extension  first  since  several  other  classes  have  been  defined  using 
similar  techniques. 

2.1_  Multivariate  I FRA.  Using  a  characterization  of  the  univariate  I FRA  class 
in  Block  and  Savits  (1976)  the  following  definition  can  be  made. 

(2.1.1)  Definition  -Let  Tm(T^ .... ,T  )  be  a  nonnegative  random  lifetime.  The 
random  vector  T  is  said  to  be  MIFRA  if 

Ea[h(T)  ]  ±E[ha(T/a)] 

for  all  continuous  nonnegative  nondecreasing  functions  h  and  all  0  <a  <1. 

This  definition  as  mentioned  above  implies  all  of  the  properties  one  would  desire 
for  a  multivariate  analog  of  the  univariate  I FRA  class.  Part  of  the  reason  for 
this  is  that  the  definition  is  equivalent  to  many  other  properties  which  are 
both  theoretically  and  intuitively  appealing.  The  statement  and  proofs  of 
these  results  are  given  below;  the  form  in  which  these  are  presented  is  influ¬ 
enced  by  the  paper  of  Marshall  and  Shaked  (1982)  who  defined  a  similar  MNBU 
class. 

Notes.  1)  Obviously  in  (2.1.1)  we  need  only  consider  h  defined  on  {x|x>0}. 

Hence  all  of  the  functions  and  sets  mentioned  below  are  assumed  to  be  Borel 
measurable  in  1R^. 

2)  We  say  a  function  g  is  homogeneous  ( subhomo geneous)  on  if 

ag(t)  -  (_<)  g(3£)  for  all  ()_<<*  <_1,  0<t. 

3)  A  is  an  upper  set  if  x«  A  and  x  implies  A. 


(2.1.2)  Theorem  The  following  conditions  are  all  equivalent  to  T  being  MIFRA. 

i)  Pa{T  e  A}  _^P{T  e  aA}  for  all  open  upper  sets  in  1r",  all  0<a^l. 

ii)  Pa{T  «  A}  _<  P{T  £  aA}  for  all  upper  sets  in  all  0<a<_l. 

(i.e.  E°($(T) )  ^E(<j>a(T/a) )  for  all  nonnegative,  binary,  nondecreasing 

b  on  JR°) . 

HI)  E^M^T) )  ^E(h(T/a) )  for  all  nonnegative,  nondecreasing  h  on  ffi", 

all  0  <  a  <c  1, 

n 

iv)  For  all  nonnegative,  nondecreasing,  subhomo geneous  h  on  ]R+, 
h(T)  is  I FRA. 

v)  For  all  nonnegative,  nondecreasing,  homogeneous  h  on  h(T)  is 

IFRA. 


Proof:  i)->ii) :  By  Theorem  3.3  of  Esary,  Proschan  and  Walkup  (1967)  for  an 

upper  set  A  and  any  c  >  0  there  is  an  open  upper  set  A^  such  that  Ac  A^  and 


P{T  £  aA^}  <_P{T  c  aA}  +e.  Thus 


p0t{T  €  A}  <_Pa{T  «  A£}<P{Te  aA£}  <_P{T  e  aA}  +  E. 

ii)  ■>  iii) :  Let  h^»  k*l,2,...  be  an  increasing  sequence  of  increasing  step 

functions  such  that  lim  h^  ■  h.  Specifically  take 

k-+« 


\<*> 


[ 


~  if  ^ih(t)  <4-  »  i“1’2 . k2k» 

2  2  2 


k  if  h(t)  >  k. 


k2 


1 

i.e.  h.  (t)  •  l  -4“  I  (t)  where  I  is  the  indicator  function  of  the  upper  set 


i-1  2k  Ai,k 


i,k 


A  ■  (tjh(t)  }.  Thus  we  need  only  prove  the  result  for  functions  of  the 

i  ,k 
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,  ti.  it.  iu  1*  «*«  IV.  at.-  it.  »L  iL  iVtU  ^  iL'tL'tl.  a 


h(,t)  -  l  a  I  (c),  a  >0, 
i-1  1  Ai  1 

where  A^,...,  Affl  are  upper  sets,  since  Che  remainder  follows  by  Che  monocone 
convergence  cheorem.  We  have 


Ea(  *  a.I.  £))-[  I  arPtT^A^]0!  [  T  {^£aV]< “ 

i*l  1  Ai  i-1 

m  f  —  m 

-  [  T  (  a“l“(c/a)dF(c)}a]ai  l  a  I  (c/ct)dF(c) 

i-1  J  1  A  "  “  i-lJ 

-  E([  I  a  I  (T/o)  ]°) 

i-1  i 

where  Che  lasc  inequaliCy  is  due  Co  Minkowski. 

iii)  ->  Def.  Obvious . 

Def.  ->  i).  From  Esary,  Proschan  and  Walkup  (1967)  for  any  open  upper  sec  A 
Chere  exist  nonnegaCive,  nondecreasing,  conCinuous  funcCions  h^  such  thaC  +  IA» 
Then  apply  che  monocone  convergence  Cheorem. 

iii)  ->  iv).  Lee  h  be  nonnegacive,  nondecreasing  and  subhomogeneous .  Then 

Pa{h(T)  >  C>  -Ea(I(tj0o)(h(T)))  lE(I(t>a>)(h(T/ct))) 

<  E  (I .  h(T) ))  -  P{  h(T)  >  a  c} 

—  (C,®)  a  —  — 

where  Che  firsc  inequaliCy  follows  from  iii)  and  che  second  by  che  subhomogeniecy . 

iv)  ->  v) :  Obvious. 

v)  ->  i).  LeC  A  be  an  open  upper  sec  and  define 


<C)  -  £ 


sup  {9  >0:  Tt_eA}  if  { 0  >  0 :  ^  teA}  i 


otherwise 


Then  h  is  nonnegacive,  nondecreasing  and  homogeneous.  Thus 
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Pa{T  £  A}  -  Pa{h(T)  >  1)  _<P(h(T)  >  a}  -P(T  €  aA}. 

(2.1.3)  Note.  The  following  two  alternate  conditions  could  also  have  been 

added  to  the  above  list  of  equivalent  conditions  (provided  F(0)  *1). 

n 

vi)  Pa{T  £  A}  <  P{T  £  aA}  for  each  set  A  of  the  fora  A  -  u  A  where 

-  -  -  i-1  1 

A.  *  {xlx  > x. },  x.  £  ]R^  and  for  all  0  <  a  <  1. 

vii)  For  each  k  »  1,2, . . . ,  for  each  a^  ,  i  ■  1, . . .  ,k,  j  ■  1, . . .  ,n,  0  ^_a^  l.00* 
and  for  each  coherent  life  function  t  of  order  kn 

t  (a^T^  ,a^2^^ » •  •  • » ain"^i  ,a21^2  *  *  *  *  *  ^kn^n^  I  FRA.  (See  Block  and 

Savits  (1980)  for  a  definition  of  coherent  life  function  and  for 
some  details  of  the  proof). 

In  conjunction  with  the  preceding  result  the  following  lemma  makes  it  easy 
to  demonstrate  that  a  host  of  different  lifetimes  are  MIFRA. 

(2.1. A)  Lemma.  Let  T  be  MIFRA  and  i h . tfc  be  any  continuous.  subhomoKeneous 

— “ 1  m 

functions  of  n  variables.  Then  if  S  -ijj.(T)  for  i«l,...,m,  S-(S,,...,S  )  is 

11  —  —  l  m 

MIFRA. 

Proof :  This  follows  easily  by  considering  a  nonnegative,  increasing,  continuous 
function  h  of  m  variables  and  applying  the  MIFRA  property  of  T  and  the  monoto¬ 
nicity  of  the  ij 

(2.1.5)  Corollary.  Let  x  ,...,x  be  coherent  life  functions  and  T  be  MIFRA.  Then 

x  in  — 

(T1(T),...,xm(T))  is  MIFRA. 

Proof:  Since  coherent  life  functions  are  homogeneous  this  follows  easily. 

(2.1.6) .  Example.  Let  X^...^  be  independent  I  FRA  lifetimes  and  4>  i  Sj^  c  (1,2, . . .  ,n} , 
i  ■  1, • • . ,m.  Since  it  is  not  hard  to  show  that  independent  I FRA  lifetimes  are 

MIFRA,  it  follows  that  T  -  min  X  i«l,...,m  are  MIFRA.  Since  many  different 

J£St  J 

types  of  multivariate  IFRA  can  be  generated  in  the  above  way,  the  example  shows 
that  any  of  these  are  MIFRA.  See  Esary  and  Marshall  (1979)  where  various  types 
of  multivariate  IFRA  of  the  type  in  this  example  are  defined.  See  Block  and 
Savits  (1982)  for  relationships  among  these  various  definitions. 


Multivariate  shock  models  with  multivariate  IFRA  properties  have  been 
treated  in  Marshall  and  Shaked  (1979)  and  in  Savits  and  Shaked  (1981). 

2.2  Multivariate  NBU.  As  with  all  of  the  multivariate  classes,  the  need  for 
each  of  them  is  evident  becuase  of  the  usefulness  of  the  corresponding  univariate 
class.  The  only  difference  is  that  in  the  multivariate  case,  the  Independence  of 
the  components  is  lacking.  In  particular  the  concept  of  NBU  is  fundamental  in 
disucssing  maintenance  policies  in  a  single  component  system.  For  a  multicomponent 
system,  where  components  are  dependent,  marginally  components  satisfy  the  univariate 
NBU  property  under  various  maintenance  protocols.  However,  a  joint  concept  de¬ 
scribing  the  interaction  of  all  the  components  is  necessary.  Hence  multivariate 
NBU  concepts  are  required. 

Mose  of  the  earliest  definition  of  multivariate  NBU  (see  for  example  Buchanan 
and  Singpurwalla  (1977))  consisted  of  various  generalizations  of  the  defining 
property  of  the  univariate  NBU  class.  For  a  survey  of  these  see  definitions  (1)- 
(5)  of  Section  5  of  Block  and  Savits  (1981).  For  shock  models  which  satisfy  these 
definitions  see  Marshall  and  Shaked,  (1979),  Griffith  (1982),  Ebrahlmi  and  Ghosh 
(1981)  and  Klefsjo  (1982),  Other  definitions  involving  generalizations  of  pro¬ 
perties  of  univariate  NBU  distributions  are  given  by  (7)-(9)  of  the  same  reference. 
These  are  similar  to  definitions  used  by  Esary  and  Marshall  (1979)  to  define  multi¬ 
variate  I FRA  distributions.  Properties  (7)  and  (8)  of  the  Block  and  Savits  (1981) 
reference  represent  a  certain  type  of  definition  and  bear  repeating  here.  The 
vector  T  is  said  to  be  multivariate  NBU  if: 

t(T, ,...,T  )  is  NBU  for  all  t  in  a  certain  class 
1  n 

of  life  functions;  (2.2.1) 

There  exist  independent  NBU  X^,...,X^  and  life  functions 
r. ,  i-l,...,n  in  a  certain  class  such  that  T  -t,(X), 
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El-Neweihi,  Pros chan  and  Sethuraman  (1983}  have  considered  a  special  case  of 
(2.2.2)  where  Che  are  minimums  and  have  related  this  case  co  some  other 
definitions  including  the  special  case  of  (2.2.1)  where  t  is  any  minimum. 

As  shown  in  Theorem  2.1, definitions  involving  increasing  functions  can  be 
given  equivalently  in  terms  of  upper  (or  open  upper)  sets.  Two  multivariate  NBU 
definitions  which  were  given  in  terms  of  upper  sets  were  those  of  El-Neweihi 
(1981)  and  Marshall  and  Shaked  (1982).  These  are  respectively: 

For  every  upper  set  Ac  1R+  and  for  every  0  <  p  <  1 

t  J** 

P(l£  A}  <  P(jnin(-’  A)  (2.2.3) 

—  —  n 

where  T,TT' ,T"  are  independent  and  have  the  same  distribution. 

For  every  upper  set  A  and  for  every  a  >  0,  8  >  0 

P{T  €  (a+8)A}  <  P(T  €  aA}  P{Te8A}.  (2.2.4) 

Relationships  among  these  definitions  are  given  in  El-Neweihi  (1981).  A  more 
restrictive  definition  than  either  of  the  above  has  been  given  in  Berg  and 
Res  ten  (1984): 

For  every  upper  A,! c  k  , 

P(T  «  A«-B)  ^  P(T<  A)P(T<  B)  (2.2.5) 

This  definition  was  snown  to  be  useful  in  percolation  theory  as  well  as  relia- 
o  i .  it  v  tneorv  . 

a  general  framework  involving  generalisations  of  the  concept  (2.2.1)  called 
taxing  :  ne  -closure  it  f  and  of  the  concept  (2.2.2)  called  C-generating  from 
*  utter*  1  .»  the  .am  of  univariate  NBU  lifetimes  in  (2.2.1)  and  (2.2.2)) 
riA*  aeer  g:  •>*  Marshall  and  Shaked  (19S<J.  Many  of  the  previous  NBU  defi¬ 

cit  ion*  are  >r gen  lied  within  this  framework.  A  similar  remark  applies  when 
the  . aaaee  f  are  exponential,  ! FR .  I FUA  and  NBUE.  See  Marshall  and  Shaked 


2.3  Multivariate  NBUE.  Along  with  the  multivariate  NBU  versions  of  Buchanan 
and  Singpurwalla  (1977)  are  integrated  versions  of  these  definitions.  These 
authors  give  three  versions  of  multivariate  NBUE.  The  relations  among  these 
and  closure  properties  are  discussed  in  Ebrahimi  and  Ghosh  (1981).  Furthermore 
the  latter  authors  relate  these  multivariate  NBUE  definitions  to  four  defi¬ 
nitions  of  multivariate  NBU  (i.e.  definitions  (l)-(4)  of  Section  5  of  Block  and 
Savits  (1981)). 

Some  other  multivariate  NBUE  classes  are  mentioned  by  Block  and  Savits 
(1981)  and  Marshall  and  Shaked  (1984).  One  extension  of  a  univariate  charac¬ 
terization  of  the  NBUE  class  mentioned  in  Block  and  Savits  (1978)  has  been 
proposed  by  Savits  (1983b). 

2.4  Multivariate  IFR.  Perhaps  the  most  important  univariate  concept  in  reliability 
is  that  of  increasing  failure  rate.  One  reason  for  this  is  chat  in  a  very  simple 
and  compelling  way  this  idea  describes  the  wearout  of  a  component.  Many  engineers, 
biologists  and  actuaries  find  this  description  fundamental.  The  monotonicity  of 
the  failure  rate  function  is  simple  and  intuitive  and  occurs  in  many  physical 
situations.  This  also  is  crucial  in  the  multicomponent  case  where  the  components 
are  dependent. 

Several  authors  have  attempted  to  describe  Che  action  of  the  failure  rates 
increasing  for  n  components  simultaneously.  These  cases  were  discussed  in  Block 
and  Savits  (1981)  and  in  the  references  contained  therein. 

A  recent  definition  of  multivariate  IFR  was  given  by  Savits  (1983a)  and  is 
in  the  spirit  of  the  classes  defined  by  Block  and  Savits  (1980)  and  Marshall  and 
Shaked  (1982).  For  shock  models  involving  multivariate  IFR  concepts  see 
Ghosh  and  Ebrahimi  (1981). 


It  is  shown  in  Savits  (1983a)  that  a  univariate  lifetime  T  is  IFR  if  and 
only  if  E[h(x,T)J  is  log  concave  in  x  for  all  functions  h(x,t)  which  are  log 
concave  in  (x, t)  and  are  nondecreasing  in  t  for  each  fixed  x  >_  0.  This  leads 
to  the  following  multivariate  definition. 

(2.4.1)  Definition.  Let  T  be  a  nonnegative  random  vector.  Then  T  has  an  MIFR 
distribution  if  E[h(x,T)]  is  log  concave  in  x  for  all  functions  h(x,t_)  which 
are  log  concave  in  (x,0  and  nondecreasing  in  t  >_0  for  each  fixed  x^O. 

This  class  enjoys  many  closure  properties.  Among  these  are  that  all  mar¬ 
ginals  are  MIFR,  conjunction  of  independent  MIFR  are  MIFR,  convolutions  of  MIFR 
are  MIFR,  scaled  MIFR  are  MIFR,  nonnegative  nondecreasing  concave  functions  of 
MIFR  are  MIFR,  and  weak  convergence  preserves  MIFR.  See  Savits  (1983a)  for 
details.  From  these  results  it  follows  that  the  multivariate  exponential  dis¬ 
tribution  of  Marshall  and  Olkin  (1967)  is.  MIFR,  as  are  all  distributions  with 

log  concave  densities.  Since  the  multivariate  folded  normal  has  a  log  concave 
density,  it  also  is  MIFR. 

The  technique  used  in  Definition  2.4.1  for  the  MIFR  class  extends  to  other 
multivariate  classes.  In  particular,  if  we  replace  log  concave  with  log 
subhomogeneous,  we  get  the  same  multivariate  I FRA  class  as  in  Definition  2.1.1; 
if  we  replace  log  concave  with  log  subadditive,  we  get  a  new  multivariate  NBU 
class  which  is  between  that  of  (2.2.3)  and  (2.2.4).  For  more  details  see 
Savits  (1983a,  1983b). 

2*5  Multivariate  DMRL  and  HNBUE.  Few  definitions  of  multivariate  DMKL  have  been 
discussed  in  the  literature,  although  E.  El-Neweihi  has  privately  communicated 
one  to  us.  Since  developments  are  premature  with  respect  to  this  class  we 
will  not  go  into  details. 
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Multivariate  extensions  of  the  HNBUE  class  have  been  proposed  by  Basu  and 
Ebrahimi  (.1981)  and  Klefsjo  (1980).  The  extensions  of  the  former  authors  are 
similar  in  spirit  to  the  multivariate  NBUE  classes  of  Ghosh  and  Ebrahimi  (1981) 
The  latter  author's  definition  extend  the  univariate  definition  by  replacing 
the  univariate  exponential  distribution  with  the  bivariate  Marshall  and  Olkin 
(.1967)  distribution  and  considering  various  multivariate  versions  of  the  defi¬ 
nition. 

Basu  and  Ebrahimi  (1981)  show  relationships  among  their  definitions  and 
Klefsjo's,  give  some  closure  properties  and  also  point  out  relations  with 
multivariate  NBUE  classes. 
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